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Background: The spatial distribution of charge and magnetization in the proton and neutron 
are encoded in the nucleon electromagnetic form factors. The form factors are all approximated 
by a simple dipole function, normalized to the charge or magnetic moment of the nucleon. The 
differences between the proton and neutron form factors and the deviation of Q from zero are 
sensitive to the difference between up- and down-quark contributions to the form factors. 

Purpose: Recent measurements of G E up to 3.4 (GeV/c) 2 allow for a much more detailed 
examination of the form factors. The flavor-separated form factors provide information on the 
quark flavor dependence of the nucleon structure and test theoretical models of the form factors. 

Methods: We combine recent measurements of the neutron form factors with updated extractions 
of the proton form factors, accounting for two-photon exchange corrections and including an estimate 
of the uncertainties for all of the form factors to obtain a complete set of measurements up to 
Q 2 « 4 (GeV/c) 2 . We use this to extract the up- and down-quark contributions which we compare 
to recent fits and calculations. 

Results: We find a large differences between the up- and down-quark contributions to Ge and 
Gm, implying significant flavor dependence in the charge and magnetization distributions. The 
rapid falloff of the ratio G F E / G V M does not appear in the individual quark form factors, but arises 
from a cancellation between the up- and down-quark contributions. We see indications that the 
down-quark contributions to the Dirac and Pauli form factors deviate from the suggested 1/Q 4 
scaling behavior suggested by a previous analysis. Finally, while recent models provide a generally 
good qualitative description of the data, the down-quark contribution to Ge/Gm and F2/F1 are 
not reproduced by any of the models. 

PACS numbers: 25.30.Bf, 13.40.Gp, 14.20.Dh 



I. INTRODUCTION 

The nucleon electromagnetic form factors provide in- 
formation on the spatial distribution of charge and mag- 
netization of the nucleon corresponding to a Fourier 
transform of the nucleons charge or magnetization den- 
sity in a non-relativistic picture. As such, the form fac- 
tors provide some of the most direct constraints on the 
partonic structure of the nucleon. The form factors are 
functions only of the four-momentum transfer squared, 
Q 2 , and can also be expressed in terms of the Dirac, 
Fi(Q 2 ), and Pauli, F^Q 2 ), form factors, which are re- 
lated to the electric and magnetic form factors: 

G e (Q 2 ) = F 1 (Q 2 )-tF 2 (Q 2 ), 
G m (Q 2 )^F 1 (Q 2 ) + F 2 (Q 2 ), (1) 

where r = Q 2 /4M^ and Mjv is the mass of the nucleon. 
In the limit Q 2 — > 0, Ge and Gm become the charge and 
magnetic moment of the nucleon, while F\ and F2 yield 
the charge and anomalous magnetic moment. 

It has long been known that G p El G P M , and approx- 
imately follow the dipole form, Go = (l+QVQo) -2 
Qo = 0.71 (GeV/c) 2 , up to Q 2 =5-10 (GeV/c) 2 , while the 
neutron electric form factor, G E , is close to zero. This 
observation is consistent with the simple, non-relativistic 
interpretation in which the charge and magnetization of 
the nucleon is carried by the quarks, and the up and the 



down quark have similar spatial distributions. This yields 
identical contributions for all form factors except for G E , 
for which there is a nearly complete cancellation between 
the up- and down-quark charge distributions. As such, 
measurements of G E played an important role, demon- 
strating that there is a measurable difference between the 
up- and down-quark distributions. This is qualitatively 
consistent with the pion cloud picture of the neutron, 
where a positive core and negative cloud arise from a 
virtual p+7r~ component of the neutron structure. 

Details of the recent progress in measurements of the 
nucleon electromagnetic form factors can be found in re- 
cent global analyses and reviews 0-0|- Several of these 
new measurements demonstrate the limitations of the 
simple, nonrelativistic picture. The decrease of the ratio 
G E jG V f, with Q 2 , as observed in polarization measure- 
ments [a-Hlj. provided a clear demonstration that the 
form factors were not simply the sum of dipole-like con- 
tributions from the up and down quarks. While these 
results were inconsistent [d3 - [14T | with earlier extractions 
based on Rosenbluth separation techniques [ll|, this is 
now widely believed to be the result of small two-photon 
exchange (TPE) corrections which yield a small angu- 
lar dependence to the cross section [l6|, [TtJ , mimicking 
the small signal expected from the contribution of G E at 
high Q 2 . More recent measurements of the neutron form 
factors [HI, EH have provided a complete set of data up 
to Q 2 =3A (GeV/c) 2 , which has enabled a detailed com- 
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parison of the up- and down-quark contribution from the 
high-Q 2 G E measurements, as well as a detailed compar- 
ison of the proton and neutron magnetic form factors. 

Assuming isospin and charge symmetry and neglecting 
the contribution of strange quarks allows us to express 
the nucleon form factors in terms of the up- and down- 
quark contributions [2(| [2l| , 
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This yields the following expression for the up- and down- 
quark contributions to the proton form factors: 
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with similar expressions for F± and F2. In this con- 
vention, G E M represents the up-quark distribution in 
the proton and the down-quark distribution in the neu- 
tron. Because the charge is factored out from the up- 
and down-quark contributions, the Q 2 = values for 
these form factors are G E = 2, G E = 1, while the 
quark magnetic moments are taken to be the Q 2 = 
limit of the magnetic form factors: \i u = (2fi p + fj, n ) and 
Md = (Mp + 2/in). Note that the up- and down-quark con- 
tributions as defined here are the combined quark and 
antiquark contributions, and so represent the difference 
between the quark and antiquark distributions, due to 
the charge weighting of the quark and antiquark contri- 
butions to the form factors. 



II. FORM FACTOR INPUT AND 
TWO-PHOTON EXCHANGE CORRECTIONS 

Recently, the ratio R n — /j, n G E /G M of the 
neutron was measured at Jefferson Lab up to 
Q 2 =3.4 (GeV/c) 2 Ql]. These data, combined FL = 
fi p G E /G p M measurements in the same Q 2 range [9l4ll1| 
allowed for the first time a comparison of the behavior of 
F% I F™ and F p /Ff, as well as a separation of the up- and 
down-quark contributions to the form factors. In the pio- 
neering work of [7| , which will be referred to as "C JRW" 
throughout this text, measurements of R n = n n G E /G M 
for 0.30 < Q 2 < 3.20 (GeV/c) 2 were combined with pa- 
rameterizations of G E , G p Ml and G M [22j to examine the 
flavor-separated contributions and the ratio 
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In the CJRW analysis, only the uncertainty from R n was 
included in the analysis, as this was the largest source of 
uncertainty in the quantities they examined. 

Refs. @, Qjl] provided the first results for flavor- 
separated form factors at high Q 2 values, and demon- 
strated a significant difference between the up- and down- 
quark contributions. We expand on their analysis mainly 



be accounting for two effects that were not included in 
their initial result. First, we include uncertainties asso- 
ciated with all of the form factors, as they are an impor- 
tant contribution for some of the extracted quantities. 
In addition, the proton form factor parameterization (22| 
used in their analysis did not apply any two-photon ex- 
change corrections, although data was selected with an 
eye towards reducing the impact of TPE corrections. We 
address this by using an extraction of the proton form 
factors which includes TPE corrections. 

We use Kelly's parameterization of G^ along with the 
fit to R n — (inGg/G 7 ^ from Riordan et al. Ta king 
the full error correlation matrix for R n and G^ [HI, [24| 
yields uncertainties on G E that are significantly smaller 
than the uncertainties on the individual measurements, 
due to the simple functional form of the R n parameteri- 
zation. To account for this, we scale up the uncertainty 
on G E by a factor of two to provide more realistic uncer- 
tainties in the flavor-separated results. 

The leading TPE effect on the electron-proton elastic 
scattering cross section, a ep , comes from the interference 
of the one- and two-photon exchange amplitudes which 
yield a small correction to both the cross section and 
recoil polarization measurements. The angular depen- 
dence of this correction to the cross section can include 
a much larger effect on the extracted form factors [25| . 
while the recoil polarization data do not have a simi- 
lar amplification of the effect. Recent measurements of 
the angular dependence of R p also suggest small TPE 
contributions to the polarization measurements [26| . To 
account for the TPE contribution to a epi one can add 
an additional term which forces the Rosenbluth extrac- 
tion to yield the same value of G P E /G P M as the polar- 
ization transfer data. We account for TPE effects by 
using the extraction of G P E and G P M from an analysis 
which constrains the TPE corrections based on the dis- 
crepancy between Rosenbluth and polarization measure- 
ments . We use the form factors extracted based on 
the TPE parametrization from Borisyuk and Kobushkin 
(BK parametrization) [28[ , which takes the corrections to 
be linear (29| in the virtual photon polarization parame- 
ter, e, and constrains the correction to vanish in the limit 
of small angle scattering, as expected from charge conju- 
gation and crossing symmetry [l7l |30| , and as observed 
in comparisons of positron-electron scattering [3l[ . 

In this analysis, we take the extraction of G P E and G P M 
from Ref. 12711 for several electron-proton scattering mea- 
surements [H, l32T - l37j . and add data at lower Q 2 from 
Ref. [38| analyzed following the same procedure. This 
provides values of G P E and G P M ,with TPE corrections con- 
strained by polarization transfer data. In the analysis 
of Ref. [23], no uncertainty is applied for the parame- 
terization of R p = G P E /G P M . For our analysis, we in- 
clude an additional fractional uncertainty on R p equal to 
5r = 0.01Q 2 (Q 2 in (GeV/c) 2 ) to account for the uncer- 
tainty in the polarization transfer ratio R. 

This approach to constraining TPE is not expected to 
be as reliable at low Q 2 values, as the difference between 
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the two techniques is significantly smaller in this region, 
and because the parameterization of R p from polarization 
data included very little low-Q 2 data. However, recent 
low-Q 2 measurements [39U42] ] suggest that the parame- 
terization is relatively reliable for the range of data exam- 
ined here, and the agreement between these polarization 
measurements and new Rosenbluth separation data (43[ 
support the idea that the corrections are relatively small. 
Thus, the final result should be somewhat insensitive to 
the exact details of the extraction of the TPE effects in 
this region. Because the TPE corrections are still impor- 
tant in the limit of low Q 2 [l?], 143 - 146} , we compare our 
results to an extraction of the flavor-separated form fac- 
tors using theproton form factor parameterization from 
Refs. [4j and [fjj. Both of these extractions include TPE 
corrections calculated in a hadronic framework [U |48j] 
which is expected to be more reliable at low values of 
Q 2 , allowing for a check of our low-Q 2 phenomenological 
TPE extraction. The fit of Venkat et al. also includes 
additional polarization data, in particular at low Q 2 val- 
ues [391—42] , and includes a more careful evaluation of the 
low-Q 2 behavior of the fit. 

Note that a more complete flavor separation at the low- 
est Q 2 values would involve an updated extraction of the 
form factors, including calculated TPE corrections and 
the most recent form factor data [39M43I ] , along with con- 
straints on strange-quark contributions taking measure- 
ments of elastic parity- violating electron scattering [49l — 
[62J. For this work, the primary focus is at somewhat 
higher Q 2 data, and the comparison of our TPE correc- 
tions to the hadronic corrections [HI applied in the 
recent proton fits [U, Q should provide an idea of the 
robustness of the low Q 2 results. 



III. RECENT THEORETICAL PREDICTIONS 

In this section we summarize several recent theoretical 
studies of the nucleon elastic form factors which we will 
compare to our extracted flavor-separated form factors. 

Cloet et al. [631 ) presented a calculation of a dressed- 
quark core contribution to the nucleon electromagnetic 
form factors defined by the solution of a Pioncare covari- 
ant Faddeev equation. This calculation includes dressed- 
quark anomalous magnetic moment within the frame- 
work of Dyson-Schwinger equations (DSEs). The Fad- 
deev equation was described by specifying that quarks 
are dressed, and two of the three dressed quarks are 
always correlated as color-3 diquarks. The nucleon is 
represented by a Faddeev amplitude of the form \P = 
("ti + ^2 + ^3) with ^1.2 are obtained from ^3 by 
a cyclic permutation with ^^(pi, on, tj) expressed as 
sum of scalar- and axial- vector-diquark correlations with 
(pi,Oii,Ti) being the momentum, spin, and isospin labels 
of the quarks. The Faddeev equation satisfied by ^3 
was constructed by specifying the dressed-quark propaga- 
tor, diquark Bethe-Salpeter amplitudes, and the diquark 
propagators. The nucleon-photon vertex was calculated 



using six diagrams, with photon coupling to the quark or 
the diquark along with loop and exchange terms. In this 
approach, the Faddeev equation has only two new pa- 
rameters in the nucleon sector: the masses of the scalar 
and axial-vector diquarks. The scalar mass is set by re- 
quiring a nucleon mass Mjq = 1.18 GeV, and the axial- 
vector mass is chosen so as Ma = 1.33 GeV. The proton 
has a mass larger than the physical value to allow for an 
additional contributions from the pseudoscalar mesons. 
The quark, diquark, and exchange (two body) contribu- 
tions to the nucleons form factors were calculated up to 
Q 2 =12 (GeV/c) 2 . In addition, the decomposition ac- 
cording to diagram diquark content and their flavor con- 
tributions were also calculated for the same Q 2 range. 
The DSE approach aims to simultaneously describe me- 
son and baryon physics |64| with only a few parameters, 
most of which are fixed to static properties of the pion. 
This calculation will be referred to as "DSE" through- 
out this text. Because the calculation does not include 
pion cloud contributions, the masses and magnetic mo- 
ments are not expected to reproduce the physical values. 
In Ref . [63| , comparisons to data were made in terms of 
y = Q 2 /M 2 . In our analysis, we evaluate their parame- 
terization in terms of y using the physical nucleon mass, 
and the experimental magnetic moments of the proton 
and neutron. Rescaling to the physics magnetic moments 
removes the discrepancies at Q 2 — 0, although to the ex- 
tent that the calculation leaves room for additional pion 
contributions, more important at low Q 2 , this would be 
expected to worsen the agreement at larger Q 2 values. 



Cloet and Miller |65[ proposed a relativistic constituent 
quark model which is constrained by the nucleon form 
factors but also reproduces the quark spin content of the 
nucleon. This is an extension of a previous light-front 
calculation which included three constituent quarks [66j 
and predicted the falloff of G V E /G V M . In this model, it 
is assumed that the quarks are moving in a cloud of 
pions. The valance quarks are represented by quark- 
diquark combination and treated in a way consistent with 
Pioncare invariance. Due to the long range nature of 
the quarks interactions as mediated by a single pion ex- 
change, a pion emitted by a nucleon can also be absorbed 
by the same nucleon, allowing for a pion cloud contribu- 
tion. The light-front wave function that describes the 
interaction of a quark and diquark to form a nucleon 
is used to construct the Fock state needed to represent 
the nucleon. The quark-diquark approximation includes 
both scalar and axial- vector correlations, and the fla- 
vor couplings were added to obtain a symmetric spin- 
flavor wavefunction. The pion component introduced 
using a single pion loop around the bare nucleon with 
three main diagrams. These diagrams include photon 
coupling to the bare nucleon and coupling to the nucleon 
or pion in the pion loop. Terms involving jN — > nN 
couplings and the effects of intermediate A were not in- 
cluded. The model is finally expressed in terms of ten 
parameters, representing the quark and diquark masses 
and contributions to the light-front wave function and 
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a parameter describing the high-momentum behavior of 
the pion-nucleon vertex function. The parameters were 
adjusted to provide the best fit to the nucleon form fac- 
tors [11 up to Q 2 = 10 (GeV/c) 2 . The model also yields 
a quark contribution to the proton spin which is found 
to be in agreement with experimental evaluations. This 
calculation will be referred to as pion-cloud relativistic 
constituent quark model, "PC-RCQM" . 

Gonzalez-Hernandez et al. [67j provided an interpreta- 
tion of the flavor dependence of the nucleon form factors 
in terms of Generalized Parton Distributions (GPDs). 
They incorporated the Regge contribution into GPDs 
that already apply diquark models by introducing a spec- 
tral distribution for the spectator diquark mass pn(M^). 
Inclusion of the Regge contributions is crucial to ob- 
tain the correct normalized structure functions. The 
model can be summarized in the expression F{X, £, t) = 
^Mx,m(^> C> t)Rp' a ' (A, C, t) where the flavor depen- 
dence of the nucleon form factors was attributed mainly 
to Reggeon exchange contributions and handbag (quark- 
diquark) contribution. The diquark contribution was 
later represented by the two functions H and E and 
the proton-quark-diquark vertex was parameterized us- 
ing a dipole type coupling with two fit parameters m q 
and M q . The final fit to DIS structure functions, nu- 
cleon form factors, and deep virtual Compton scattering 
were performed, one without the the new CJRW extrac- 
tions [7[, and one including these results to produce im- 
proved constraints on the flavor dependence of the GPDs. 
The flavor separated nucleon form factors, separated into 
Regge and diquark contributions, were calculated up to 
Q 2 —5 (GeV/c) 2 . For F-J", the diquark contribution dom- 
inates the Regge contribution at low Q 2 and both con- 
tributions become comparable at high Q 2 . On the other 
hand, the diquark and Regge contributions are compa- 
rable for Ff and the Regge contribution dominates F% 
and F2, in particular at high Q 2 . This calculation will 
be referred to as "GPD" throughout this text. 

Recently, Rohrmoser et al. (68| analyzed the flavor 
decomposition of the nucleon electromagnetic form fac- 
tors within the framework of a relativistic constituent 
quark model whose hyperfine interaction is derived from 
Goldstone-boson exchange [69] as a result of sponta- 
neous breaking of chiral symmetry in low-energy QCD. 
In this model, nucleons are represented by three-quark- 
configuration driven by an interaction Lagrangian formed 
by coupling of Goldstone bosons with valance-quark field. 
The nucleon wavefunction has no diquark configuration 
or mesonic effects and contains non-vanishing orbital 
angular momenta and mixed-symmetric spatial wave- 
function component with relatively small non-symmetric 
contribution. The key ingredient of the nucleon wave 
functions is the interaction of the mass operator, which 
has a linear confinement, with the QCD string tension 
and chiral symmetry breaking hyperfine interaction. We 
show here the latest theoretical results [7(| and not the 
originally published [68[ and note that the two results 
are quite different. The form factors and their flavor 



decomposition were calculated up to Q 2 =4 (GeV/c) 2 . 
Throughout this text this calculation will be referred to 
as "GBE-RCQM" . 

Note that the GBE-RCQM and DSE calculations arc 
predictions for the flavor-separated form factors, as they 
do not adjust parameters to match the proton or neu- 
tron form factors, while the PC-RCQM and GPD curves 
are fits to the data using parameterizations based on the 
model. The PC-RCQM model has a total of 10 param- 
eters, while the GPD fit has 16 parameters for each of 
the GPDs, which are fit to reproduce both the form fac- 
tors and parton distributions, with additional constraints 
from lattice QCD. Thus, one expects these models to bet- 
ter reproduce the data, although this does not represent 
as conclusive of a test as it does for the DSE and GBE- 
RCQM predictions. The models which are adjusted to 
reproduce the form factor data must also be compared to 
other measurements, as mentioned above and discussed 
in more detail in the original works [65l . |67[ . 



IV. RESULTS AND DISCUSSION 

In this section, we present the flavor-separated results 
for the proton and neutron form factors. We compare the 
results to the CJRW extraction and to extractions based 
on the recent proton parameterizations from Venkat et 
al. @ ("VAMZ"), and from Arrington et al. ("AMT"), 
combined with the neutron form factors parameteriza- 
tions from Kelly [22| and Riordan et al. [T3. Note that 
these are not independent extractions; they are all based 
on fits to relatively up-to-date sets of form factor mea- 
surements. The comparison to the CJRW result allows 
us to examine the impact of the TPE corrections to the 
proton form factors, as well as the uncertainties in the 
extraction of G P E , G P M , and G 1 ^ , which were neglected in 
the CJRW analysis. The comparisons to the AMT and 
VAMZ fits provide sensitivity to the TPE corrections at 
low Q 2 , as discussed in the previous section, as well as 
the impact of recent polarization transfer data [4l|, |42j , 
which are only included in the VAMZ result. In addition, 
we compare the results to a set of recent nucleon form 
factor models, described in Sec. IHII 

Figures [T] and [5] show the proton Sachs form factors 
and their contributions from up and down quarks. Our 
extracted values are included in the online supplemen- 
tal material [7l[. The top panels show the proton form 
factors used in the extraction, normalized to the dipole 
form, along with the values from the CJRW analysis, the 
AMT and VAMZ fits, and recent calculations. 

Note that for G P E and G P M , as well as the flavor- 
separated Gm values, the CJRW points have no uncer- 
tainty. This is simply because they extract these directly 
from parameterizations of G E , G P M , and G^ and do not 
include any uncertainty in the fits. Obviously, when look- 
ing at form factors that are independent of the R n mea- 
surements, these uncertainties cannot be neglected. 

Our results are otherwise in relatively good agreement 
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FIG. 1: (color online) G p e /Gd [top] and its flavor-separated 
contributions 2/3(G|/G D ) [middle;}, and -1/3(G%/G D ) [bot- 
tom] as obtained from Refs. [l^ . [32433] based on our fit and 
the CJRW extractions 0]. Also shown are the AMT [3 and 
VAMZ fits % and the values from the GBE-RCQM [gg,^. 
PC RCQM (i.Y. the DSE :63], and the GPD [67[ models. 



with the CJRW analysis. The different treatment of 
TPE corrections in the proton form factor input yields a 
small difference in G P E and G P M for Q 2 values near 0.5- 
1.5 (GeV/c) 2 . The CJRW result are in better agreement 
with the improved treatment of TPE corrections in the 
AMT and VAMZ fits for G E , while our extraction is in 
better agreement for G P AI . At larger Q 2 values, both anal- 
yses and the global fits yield consistent results. 

As expected, the up-quark contribution dominates 
both the charge and magnetic form factor for the pro- 
ton. Examining the contributions to G P E , we see that 
both the up- and down-quark contributions have signifi- 
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FIG. 2: (color online) G p m /Gd and its flavor-separated con- 
tributions. Curves and data points are the same as in Fig. [T] 



cant deviations from the dipole form. At low Q 2 values, 
the increase in G e /Gd is compensated by a decrease in 
G e /Gd, yielding a small Q 2 dependence in G p e /Gd- At 
higher Q 2 values, G d E is consistent with the dipole form, 
and the decrease of G e /Gd leads to the overall falloff in 
G p e /Gd- Note that the cancellation between the positive 
but slowly decreasing value of G E / Gd and the negative 
but nearly constant value of G e /Gd enhances the Q 2 de- 
pendence seen in the up-quark contribution. Thus, the 
rapid linear falloff observed in G P E /G P M is a result of the 
cancellation between the contributions from the up and 
down quarks. This behavior is therefore connected to the 
difference in the up and down quark distributions, rather 
than the overall shape of the quark distributions. 

For G P M , both up and down quarks have smaller devi- 
ations from the dipole form. At very low Q 2 values, both 
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FIG. 3: (color online) G^/Gd and its flavor-separated con- 
tributions. Note that G% (middle figure) represents the up- 
quark contribution to the proton, and therefore the down- 
quark contribution in the neutron, and so is multiplied by the 
charge of the down quark. Curves and data points are the 
same as in Fig. \T\ 
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FIG. 4: (color online) G^m/Gd and its flavor-separated con- 
tributions. Note that G\j (middle figure) represents the up- 
quark contribution to the proton, and therefore the down- 
quark contribution in the neutron, and so is multiplied by the 
charge of the down quark. Curves and data points are the 
same as in Fig. [T] 



have a slightly increasing contribution, yielding a roughly 
10% increase in G P M between Q 2 of and 1 (GeV/c) 2 . 
Above this, the slow increase in and slow decrease 
in Gfj yielding a near-perfect agreement of G V M with the 
dipole form up to Q 2 —4.5 (GeV/c) 2 , even though both 
the up and down contributions have significant devia- 
tions. 

Figures [3] and 2] show the Sachs form factors of the 
neutron along with their breakdown into up- and down- 
quark contributions. In this case, our results for and 
G^ simply reflect the values and uncertainties of the 
form factor parameterizations of G\ [l9[ and GJf [22j . 



except for the factor of two enhancement on SG 7 ^ we ap- 
ply to better match the experimental uncertainties. For 
the CJRW results, and the flavor-separated results 
include the uncertainties associated with the direct mea- 
surements of R n , while the values for G^j and it's up- 
and down-quark contributions are based entirely on the 
parameterizations of the proton and neutron magnetic 
form factors, with no uncertainties included. 

Unlike in the case of the proton, the up and down 
quarks both yield large contributions to the neutron form 
factors. The strong Q 2 dependence of G^/Gd at low 
Q 2 , where G%/Gd is relatively flat, yields the rise in 
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GJ, while at larger Q 2 values, G e /Gd stops rising and 
G E grows slowly compared to the dipole form due to 
the small Q 2 dependence in G%. As with the proton, 
the contributions to have small deviations from the 
dipole, although the contribution from the down quark, 
which has larger deviations, yields a small Q 2 dependence 
in Gm/Gd at larger Q 2 values. 

Figures [TpE] also show the flavor-separated contribu- 
tions from form factor parameterizations (AMT [4[ and 
VAMZ L 6i]) and the calculations discussed in Sec. IIIII As 
mentioned above, the parameterizations are fits that in- 
clude much of the data included in these extractions, and 
so at high Q 2 yield consistent results with the data. At 
low Q 2 , they help show the impact of two-photon ex- 
change corrections which are neglected in the CJRW ex- 
traction and treated in a way that is less reliable at low 
Q 2 in our analysis. 

The calculations all give a reasonable qualitative de- 
scription of the up- and down-quark contributions, show- 
ing G e /Gd rising at high Q 2 , G e /Gd rising and then 
leveling off or falling, and relatively little Q 2 dependence 
in the up- and down-contributions to Gm /Gd- The GPD 
model [67| gives the best description of the data, with 
only small deviations at large Q 2 . This is not surprising 
as it includes an essentially complete set of data in fit- 
ting the GPDs, and uses a GPD parameterizations with 
sufficient flexibility to reproduce of the data. The PC- 
RCQM result [65| also does a good job in reproducing the 
behavior of the data, although with significantly larger 
deviations in G d M (and G M ) than in the other form fac- 
tors. The GBE-RCQM calculation [H| does a fairly good 
job of reproducing the data at small low Q 2 values, but 
above 1 (GeV/c) 2 , shows large deviations from the data 
in both the flavor-separated and the proton and neutron 
form factors. While it does not reproduce the data as well 
as the GPD or PC-RCQM curves, it is a parameter-free 
calculation, making the overall agreement rather remark- 
able. 

The DSE calculation [63| has significant deviations at 
both low and high Q 2 . However, in the DSE approach, 
the nucleon mass and quark magnetic moments are not 
forced to reproduce the physical values, as is required be- 
cause additional contributions from pseudoscalar mesons, 
the "pion cloud" contributions, excluded in this calcula- 
tion, will bring these closer to the physical values. The 
missing pseudoscalar meson contributions arc expected 
to modify the behavior at lower Q 2 values, while the 
unphysical nucleon mass in this model may modify the 
comparison at higher Q 2 values. For this comparison, we 
used the physical nucleon masses and magnetic moments 
to partially account for this difference. However, while 
taking the physical magnetic moments yields the correct 
limit as Q 2 — > 0, it may worsen the agreement at larger 
Q 2 values, where the pion cloud contributions are not ex- 
pected to be as important. Note that the DSE calculation 
has most of its parameters fixed based on the calculation 
of light mesons, and the only additional parameters for 
the nucleon calculation are the diquark radii, taken to be 



commensurate with the pions charge radius. So as with 
the GBE-RCQM calculation, the result is not adjusted 
to improved agreement with the form factor data. 
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FIG. 5: (color online) The ratio Ge/Gm for the up- and 
down-quark contributions, normalized to unity at Q 2 = 0. 

Figure [5] shows the ratio Ge/^Gm for both the up 
and down quarks. As before, the extractions disagree 
somewhat at low Q 2 , with the fits including a more com- 
plete TPE treatment agreeing better with CJRW below 
1 (GeV/c) 2 , and our results above. For the up quark, 
Ge/Gm has a roughly linear falloff at large Q 2 , but 
the decrease is slower than seen for the proton. For the 
down quark, a completely different behavior is seen, with 
Ge/Gm increasing in magnitude with Q 2 . This behav- 
ior is not present in any of the calculations except for the 
GPD model, which has a sufficiently flexible parameter- 
ization of the GPDs to reproduce all of the form factors. 

We now examine the Dirac and Pauli form factors, 
with an emphasis on the high-Q 2 behavior of the flavor- 
separated contributions. While perturbative-QCD be- 
havior should set in at large enough Q 2 values, the data 
do not extend into the region where one expects these 
asymptotic predictions to be valid. Nonetheless, it has 
been observed that approximate scaling of the form fac- 
tors often sets in at lower Q 2 values. By taking out the 
predicted high-Q 2 behavior, we can more easily see differ- 
ences in the Q 2 dependences of the various contributions 
to the form factors. 

Figure [5] shows the proton Dirac and Pauli form fac- 
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FIG. 6: (color online) Q 4 F[ [top], Q 6 F 2 P [middle], and 
Q 2 Fl Y-FT [bottom] from polarization measurements of R p [g, 
OP . [ill ] and G parameterization of Kelly [22] • Also shown 
are the AMT 0] and VAMZ 0] fits, and the calculations dis- 
cussed in Sec . IIHI We also show the modified pQCD scaling 
fit from Ref. [zl] with A=300 MeV, labeled "BJY-pQCD". 
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FIG. 7: (color online) Q 4 F? [top], Q 6 F? [middle], and 
Q 2 F?/F? [bottom] from polarization measurements of 
Rn [la . [73477] ] and the G 7 ^ parameterization of Kelly [2^ ]. 
Also shown are the Kelly fit [22] and the calculations pre- 
sented in Sec. IIHI 



tors and their ratio. In all cases, the leading perturba- 
tive ("pQCD") Q 2 dependence [78| is removed by scal- 
ing the results by powers of Q 2 . While pQCD suggests 
scaling behavior of the form Ff oc Q~ 4 , -Ff oc Q~ 6 , and 
Ff /Ff oc Q~ 2 , the data clearly do not support such scal- 
ing as both Q 4 Ff and Q 6 Ff increase with Q 2 . While 
Rosenbluth extractions which did not include TPE cor- 
rections observed scaling behavior in the flattening of the 
ratio Q 2 F% [F f ', the high-Q 2 recoil polarization measure- 
ments [§L l79| show that the pQCD scaling behavior is 
not observed. Note that while both Ff and Ff deviate 
from the scaling predictions, these deviations are differ- 



ent enough that the ratio also deviates from the pQCD 
expectation. 

Figure [7] shows the same quantities for the neu- 
tron. Again, neither F± nor F2 are consistent with the 
pQCD predictions. However, unlike the proton, the 
ratio F2/F1 is consistent with a constant value above 
Q 2 =1.5 (GeV/c) 2 , as noted in the CJRW analysis Q, 
although the precision of the data do not set tight con- 
straints on the Q 2 dependence. 

The DSE and GBE-RCQM calculations generally show 
similar deviations from experiment for the Dirac and 
Pauli form factors (for both the proton and the neutron), 
yielding better agreement with the ratio than the indi- 
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FIG. 8: (color online) The up-quark [top] and down-quark 
[middle] contribution to the Dirac form factor multiplied by 
Q 4 , along with their ratio [bottom]. 
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FIG. 9: (color online) The up-quark [top] and down-quark 
[middle] contribution to the Pauli form factor multiplied by 
k^jQ 4 , along with the ratio F^fF? [bottom]. 



vidual form factors. This also helps explain why these 
calculations are in somewhat better agreement with Ge 
then Gm, as Ge relates to the difference between F\ and 
F2 , yielding a partial cancellation of the calculations de- 
viations from the data. At large Q 2 values, the deviations 
in the DSE calculation grow and are of the opposite sign 
for Fx and F%. However, the calculation does not adjust 
the diquark radius to better reproduce the proton form 
factor, and the deviations from the data can be signif- 
icantly improved by increasing the diquark radius [63| . 
At large Q 2 values, the behavior of the ratio is also ex- 
tremely sensitive to the Dressed-quark mass function, as 
discussed in Sec. Ill of Ref. [8(|, due to the significant 
cancellation between F\ and F2 contributions. 

We now turn to the flavor-separated contributions of 



Fi and F2. The extracted form factors are shown in 
Figs. [8] and [9j and the values are provided in the sup- 
plemental online material [7l|. Our values for F% are 
somewhat higher than the CJRW extractions at low Q 2 , 
and the impact of this difference is seen clearly in the 
F2/F2 ratio. Again, while the G E uncertainty tends to 
dominate, inclusion of the proton and uncertainties 
has a noticeable impact on some of the flavor-separated 
results. 

In the CJRW analysis [7] it was reported that 
both Ff and Frf strikingly exhibit 1/Q 4 scaling above 
Q 2 =1.0 (GeV/c) 2 , in contrast to the up-quark form fac- 
tors which continued to rise relative to the down-quark 
values. This was in agreement with the predictions for 
the moments of the generalized parton distributions re- 
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ported in Ref. 81], although these again are based on 
fits to data sets which include the nucleon form factors, 
except for the most recent G% data. Both our results 
and the global parameterizations suggest that the down- 
quark contributions are falling slightly with respect to 
the 1/Q 4 behavior, and this falloff appears to be fairly 
clear for F 2 d . However, this behavior is sensitive to the 
parameterization of Gg at Q 2 >2 (GeV/c) 2 , which are 
constrained directly only by the two data points from 
Ref. [l|| . Similarly, all of the calculations have the down- 
quark form factors falling somewhat faster than l/Q 4 , 
suggesting that the apparent scaling behavior may not 
continue to higher Q 2 . However, it is still clear that the 
down-quark contributions fall significantly more rapidly 
than the up-quark contributions at large Q 2 , yielding a 
decrease in the magnitude of F 2 / F 2 at high Q 2 . The very 
different Q 2 dependence for the u and d quark contribu- 
tions suggests that approximate 1/Q 2 scaling of F2/F1 
for the neutron is only approximate, and may be unre- 
lated to the predicted scaling behaviors. 

The faster falloff of the down-quark contributions was 
interpreted in Ref. i7| and references therein as an in- 
dication of the possibility of sizable nonzero strange 
matrix elements at large Q 2 or the importance of di- 
quark degrees of freedom. While existing measure- 
ments of parity-violating elastic scattering yield very 
small contributions from the strange quarks up to Q 2 « 
1 (GeV/c) 2 [II 11, [II, 111, they still leave open the pos- 
sibility for significant contributions from G S E and G» 
which cancel in the parity- violating observables [13, [61| , 
although there are also results from Lattice QCD that the 
strange-quark contribution is small for both the charge 
and magnetic form factors 59, 82} . In the diquark model, 
the singly-occuring down quark in the proton is more 
likely to be associated with an axial-vector diquark than 
a scalar diquark, and the contributions of the axial- vector 
diquark yields a more rapid falloff of the form factor. 
The up quarks are generally associated with the more 
tightly bound scalar diquarks, yielding a harder form fac- 
tor [lllllll. 

The flavor dependence of the nucleon form factors as 
obtained in the CJRW extractions was reproduced quite 
well by incorporating the Regge contribution into Gen- 
eralized Parton Distributions (GPDs) that already apply 
diquark models [67|]- Inclusion of the Regge contribu- 
tions is crucial to obtain the correct normalized struc- 
ture functions. Therefore, the flavor dependence was at- 
tributed mainly to reggeon exchanges and quark-diquark 
contributions. For the Dirac Q 4 F™ form factor, the di- 
quark contribution dominates the Regge contribution at 
low Q 2 and both contributions become comparable at 
high Q 2 . On the other hand, for Q 4 F?, both the di- 
quark and Regge contributions are rather comparable. 
For the Pauli Q A F 2 d ^ form factors, the Regge contri- 
bution dominates that of the diquark contribution and, 
in particular, at high Q 2 . This again shows the impor- 
tance of the diquark contribution at high Q 2 , although 
in this framework, the Regge contributions are important 



in acheiving a better result at low Q 2 . These data will 
also allow a flavor separation at higher Q 2 values, where 
the pion cloud contributions, neglected or included in a 
less detailed fashion, are expected to be smaller. This 
will allow for a more direct test of the calculations of the 
three-quark core, with reduced uncertainties associated 
with the more poorly understood pion cloud contribu- 
tions. 
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FIG. 10: (color online) The ratios K^Fg/F? [top], and 
K^Fi/F? [bottom]. 



The ratios K^Fg/F? and K^Fg/F? are shown in 
Fig. [TO] The ratio n~ lF 2/ F i falls rapidly at low Q 2 , but 
the decrease is signficantly slower above Q 2 = 1 (GeV/c) 2 
(note the offset zero in the top panel of Fig. [T0|) . Our 
values are somewhat larger than those obtained in the 
CJRW extractions for Q 2 < 1.5 (GeV/c) 2 , due to the 
difference in the F 2 values. For the ratio K~ A F 2 /Ff, 
the extractions from the data and the fits to the form 
factors yield consistent results, with a slight increase in 
the ratio at low Q 2 and a nearly constant value above 
1 (GeV/c) 2 . However, all of the calculations show a very 
different behavior, with Ff falling more rapidly than F 2 d 
at large Q 2 values, leading to a rapid rise in the ratio. 
This is clearly an area where the models should be exam- 
ined more carefully, although the present measurements 
of G B do not rule out a significant rise in the ratio above 
1.5 GeV 2 . Measurements of G E planned for higher Q 2 af- 
ter the Jefferson Lab 12 GeV upgrade [83( will be critical 
in pinning down the behavior of this ratio. 
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V. CONCLUSIONS 

We have extracted the flavor-separated contributions 
to the elastic nucleon electromagnetic form factors based 
on parameterizations of the neutron form factors and 
their uncertainties, and proton form factor extractions 
that include phenomenological TPE corrections. The ex- 
traction is similar to that of the recent C JRW 0] analysis, 
but with an explicit treatment of two-photon exchange 
effects and the uncertainties on the proton form factor 
and the neutron magnetic form factors. The treatment 
of the TPE contributions yields differences in some of the 
results at low Q 2 values, up to ml. 5 (GeV/c) 2 . The addi- 
tional uncertainties generally have a small contribution, 
but are more important for some of the flavor-separated 
form factors. The extracted flavor-separated form factors 
are qualitatively reproduced by a range of models, with 
better quantitative agreement for those models which 
constain parameters by directly fitting to nucleon form 
factor data. 

The strong linear falloff with Q 2 observed in G V E / G P M 
is not present in either the up- or down-quark contribu- 
tions, but arises due to a cancellation between a weaker 
Q 2 dependence for the up quark and a negative but rela- 
tively Q 2 -independent contribution from the down quark. 
This indicates that the rapid falloff and the zero cross- 
ing expected near Q 2 — 10 (GeV/c) 2 is associated more 
with the difference between the up and down quark dis- 
tributions than by the isospin-averaged spatial density 
distributions. 

As noted in the previous analysis, the F± and F^ form 



factors show a different Q 2 dependence for up- and down- 
quark contributions, which can be shown to be a conse- 
quence of diquark degrees of freedom in several of the 
calculations. We see indications that the down-quark 
contributions to the Dirac and Pauli form factors deviate 
from the l/Q 4 scaling suggested in the CJRW analysis, 
and observe that there are small differences between the 
Q 2 dependence in F\ and Fi for both the up and down 
quark contributions. 

Finally, the up and down quarks yield very different 
contributions to Ge/Gm (and F2/F1), with the ratio de- 
creasing slowly with Q 2 for the up quark and increasing 
rapidly for the down quark. The down-quark contribu- 
tion to this ratio is not well reproduced by any of the 
calculations, and even the qualitative behavior is not re- 
produced in most approaches. Data at higher Q 2 will 
better constrain the behavior of these ratios, and allow 
for a more detailed evalation of the nucleon models in a 
region where pion cloud contributions are expected to be 
less significant. 
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